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Abstract 

We analyze collisional decoherence of atoms or molecules prepared in a coherent superposition 
of nondegenerate internal states at ultralow temperatures and placed in an ultracold buffer gas. 
Our analysis is applicable for an arbitrary bath particle/tracer particle mass ratio. Both elastic 
and inelastic collisions contribute to decoherence. We obtain an expression relating the observ- 
able decoherence rate to pairwise scattering properties, specifically the low-temperature scattering 
amplitudes. We consider the dependence on the bath particle/tracer particle mass ratio for the 
case of light bath and heavy tracer particles. The expressions obtained may be useful in low- 
temperature applications where accurate estimates of decoherence rates are needed. The results 
suggest a method for determining the scattering lengths of atoms and molecules in different internal 
states by measuring decoherence-induced damping of coherent oscillations. 
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I. INTRODUCTION 



Atomic and molecular interferometry experiments , precision measurements of fun- 
damental constants ^3 ^6 1 and coherent control of molecular dynamics [7] are based on gaseous 
ensembles of atoms or molecules prepared in coherent superpositions of internal (electronic, 
ro-vibrational, hyperfine, Stark or Zeeman) energy states. Atoms and molecules in coherent 
superpositions of internal states have also been proposed as building blocks for quantum 
computation and quantum information processing [8, 9]. Coherent superpositions may be 
destroyed by external field fluctuations and collisions of gas particles. Collision-induced de- 
coherence is a major limiting factor in the experimental realization of quantum computation 
and coherent control of molecular dynamics. Recent progress in the development of experi- 
mental techniques for cooling atoms and molecules to extremely low temperatures suggests 
new possibilities for precision spectroscopy measurements, coherent control of molecular 



processes and quantum computation lOj. For example, cooling molecules to ultralow tem- 
peratures allows for high-resolution spectroscopy with long interrogation times and high 
degree of control over intra- and inter-molecular interactions The translational energy 
of ultracold molecules is insignificant and can be disentangled from internal states, which can 
be exploited to develop new schemes for coherent control [12]. However, elastic and inelastic 
collisions of atoms and molecules at ultralow temperatures may be very efficient, leading to 
significant decoherence rates jl^|. In order to assess the feasibility of quantum information 
processing, quantum interferometry measurements and coherent control schemes based on 
ultracold atom and molecules, it is necessary to develop a microsopic theory of collisional 
decoherence of internal state superpositions at ultralow temperatures. 

A master equation describing decoherence of translational and internal states for a tracer 
molecule in an inert gas has been derived heuristically by Hornberger and Vacchini {l4|, and 
the Brownian motion limit of an infinitely massive tracer particle was considered by Horn- 



berger 15]. This latter result was applied to decoherence of enantiomeric states of optically 
active molecules to explain Hund's paradox, i.e. the well-known observation that chiral 
molecules are not naturally found in their achiral ground state jl6|. Although coherence in 
the tracer molecule's state in general may involve both translational and internal degrees of 
freedom, in many experiments a substantial simplification arises because the translational 
degrees of freedom are fully decohered. Internal state decoherence in this simplified case 



was studied by Vacchini 17j. Whether or not the translational degrees of freedom are fully 
decohered, coherences between internal states typically depend on the translational degrees 
of freedom. For an experiment in which only the internal state is probed, the appropriate 
statistical operator is a reduced density matrix acting on the internal state space, which 
is obtained by tracing over the translational degrees of freedom. As a consequence, the 
evolution of the internal state is typically non-Markovian and the time dependence of the 
internal state coherence is generally non-exponential jl7fl. 

The theory of collisional decoherence for a tracer particle with internal states is an ex- 



tension of the analysis of positional decoherence 
This problem has been studied in detail 
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or a tracer particle without internal states. 
-23], particularly in the Brownian motion 



limit of an infinitely massive tracer particle. Hornberger and Sipe presented a solution 
for the Brownian motion limit using a convex decomposition of the bath gas density ma- 
trix into localized wavepackets to avoid mathematical problems arising because momentum 



eigenvectors are not normalizable 



2l| . This rigorous method can be used to develop formal 



replacement rules for consistently handling the singular quantities in the momentum basis. 
The replacement rule method was extended to the case of a tracer particle with finite mass 



22|. Adler connected the replacement rules to the method used to resolve the squared Dirac 
delta function appearing in the Golden Rule derivation J^ . 

Other work on collisional decoherence of internal state superpositions has included inves- 
tigation of decoherence suppression using optical pulses {25 1 and approaches based on Monte 



Carlo simulations 
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27j. An analysis by Reinhold and co-workers predicted that measure- 



ments of decoherence-induced damping of quantum beats could be used as a sensitive probe 
of collision cross-sections J27 1. Ramakrishna and Seideman investigated decoherence of ro- 



tational wavepackets using an approach applicable for dense media 28]. This method treats 
the translational degrees of freedom phenomenologically. They also found that damping 
of coherent oscillations and relaxation in the molecular alignment may provide information 
about elastic and inelastic collision properties. 

In this paper we consider decoherence of internal state superpositions for a tracer molecule 
in the presence of an ultracold buffer gas in a trap. This process can be studied in buffer 
gas cooling experiments - an important class of experiments for the creation of ultracold 

r in sym- 



molecules and precision spectroscopy 



pathetic cooling experiments using magneto-optical traps 
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33] . The separation between 



internal energy levels of molecules is typically large compared to ultracold trap depths f^ . 
As a consequence, inelastic collisions, in which the internal state of the molecule is changed, 
release enough energy from internal states into translational motion that both the tracer 
molecule and the buffer gas atoms involved in the collision are ejected from the trap. The 
experiments begin with the translational degrees of freedom of the buffer gas and tracer 
molecule fully decohered and in thermal equilibrium at the buffer gas temperature |42|. It 
will be shown that the translational degrees of freedom remain fully decohered throughout 
the experiment. At t = the tracer molecules are prepared in a coherent superposition 
of internal states. We obtain expressions for the coherence present between internal states 
over time during the experiment in the limit of ultracold temperatures and in the limit of a 
small but finite buffer gas/tracer molecule mass ratio. These expressions give the temper- 
ature dependence of the decoherence rate at ultracold temperatures in terms of the s-wave 
scattering parameters for collisions between the buffer gas particles and the tracer molecules 
in the different internal states. 

The number of trapped molecules in experiments with ultracold gases is continually de- 
creasing due to inelastic trap loss. There are different measures of coherence between internal 
states that may be of interest depending on the experiment. In quantum computation, for 
example, the total strength of the coherent signal relative to the size of the initial ensemble 
may be important. In contrast, for an experiment which is observing coherent oscillations 
arising from interference between internal states, the quantity of interest may be coherence 
within the trapped collection of molecules, and the decreasing size of the trapped population 
may be irrelevant. We derive two different measures of internal state coherence appropri- 
ate to each of these experimental scenarios and find the temperature dependence of the 
decoherence rate at ultracold temperatures. Our analysis is based on the replacement rule 



approach of Hornberger, Sipe and Adler 21J,|24|. related to the work of Hornberger in 
Ref. [17| which provides a general theory of collisional decoherence of internal states with 
translational degrees of freedom fully decohered. Our analysis differs from Ref. [17] in that 
we incorporate the presence of inelastic trap loss to model experiments in ultracold traps, 
and we consider a particular limit of low temperatures. Our results suggest a method for 
determining scattering lengths of molecules in different states based on low-temperature 
measurements of the temperature dependence of decoherence-induced damping of coherent 
oscillations. 
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II. COLLISIONAL DECOHERENCE IN THE MOMENTUM AND INTERNAL 
STATE REPRESENTATION 



We consider an ensemble of atoms or molecules of mass M prepared at time t = 
in an internal state characterized by a reduced density matrix with elements p uu >(0) and 
placed in a bath of ultracold atoms of mass m. The quantum number v denotes an internal 
energy eigenstate with energy e v . The initial internal state may be pure or mixed. If there 
is initially coherence present between internal energy eigenstates, decoherence occurs as a 
result of collisions with buffer gas atoms. The goal of this work is to calculate the time 
evolution of the coherence between internal states v and v' using the measures of coherence 
defined in Sec. Ill CI 

We make the following assumptions to model the conditions of experiments with ultracold 

gases: 

1) At t — the bath particles and the tracer particle translational degrees of freedom are 
in thermal equilibrium at the ultracold temperature T. 

2) The tracer particle internal states \u) are nondegenerate. 

3) The buffer gas is sufficiently dilute that the time between collisions is long compared 
to the duration of a collision. This implies that we can neglect collisions involving three or 
more particles. 

4) The energy difference between the internal energy levels of the tracer particle is large 
compared to the trap depth. As a consequence, inelastic collisions transfer enough internal 
energy to translational motion to eject both the tracer particle and the buffer gas particle 
involved in the collision from the trap. 

5) The collisions between buffer gas particles keep them in thermal equilibrium. In 
addition, the previous assumption implies that there is no heating of the trapped buffer gas 
due to inelastic conditions. As a result, we may treat the buffer gas particle in every collision 
as being drawn from the thermal equilibrium ensemble at temperature T which does not 
change in time. 

6) The trapping potential is ignored: the trapped gas is treated as a uniform gas in free 
space. This is a good approximation since trap dimensions are large compared to relevant 
length scales for collisions and gas properties typically vary slowly in space. 

The tracer particle state at time t > is represented by the reduced density operator 
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pit) with matrix elements p ul/ '(P,P';t) = (Pi/|p(i)|PV) in the momentum and internal 
state basis. The bath ensemble is described by the thermal equilibrium density matrix 
corresponding to the bath temperature T, 



n 

where 



J d 3 P IpMpXpI , (i) 



g 2mk B T 

" (p) = (2*mk B T)W ' (2) 

and fl is the system box volume. Q is not to be confused with the trap volume. It is imagined 
that the free space system is placed in a box of finite volume Q and the continuum limit is 
obtained as Q — > oo. Q appears as a necessary normalization factor. 
At t — 0, the tracer particle density operator is 

m = Y, p-(o) H ¥\ « St 1 7 rflp (wFl p > < p l • < 3 > 

where p uu '{0) are the elements of the reduced density operator which describes the initial 
internal state, hence p vv > (0) = Pl>M and < p^(0)pM^)/Vp^Mp^W) < 1- 

The internal state space of the tracer particle is T-L u . The state spaces for the relative 
translational motion of the colliding pair and the motion of the center of mass are respectively 
U re i and "Hcm- 



A. Collisions without trap loss 

The Hamiltonian for the two-particle system consisting of a tracer particle and buffer gas 
particle is 

H = H + V, (4) 

where 

H = H™ + HT (5) 

is the sum of the Hamiltonians of the free particles, if™ is the Hamiltonian of the tracer 
molecule, H^ as is the Hamiltonian of the buffer gas atom and V represents the interaction 
between them. 
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At times long before and long after the collision the particles are well separated and the 
system's evolution is governed by the Hamiltonian Hq. If the colliding particles are in a pure 
quantum state, the two particle state — > e~ tHot ^ h \ ^m) as t — > — oo for some |^i n ), and 
as t — >■ oo, ^(i)) — > e~ %Hot l h l^out)- The incoming and outgoin g as ymptotic trajectories are 



related by the two-particle scattering operator S according to 



35] 



|*o»t)=5'|* 1 n) • (6) 

The duration of a collision is short compared to the time scales of the experiments. We 
may thus regard a collision occurring at time t as an effectively instantaneous transition 
between the asymptotic trajectories e ~ iHot/h |^ in ) -)• e' iH ° t/h |\EVt) = e~ iH ° t/h S \^ in ) = 
g e -iH t/h |v[/ in ^ The second equality is valid because the S operator commutes with Hq 

H. 

If the incoming two-particle state is a mixed state described by a density operator p pair , 
then the state after the collision is given by the density operator 

p pair ' = 5p pair 5 f . (7) 

The reduced one particle density operator for the tracer particle state after the collision is 
obtained by tracing over the degrees of freedom of the bath particle, 

p' = tr gas {Sr h S^} . (8) 

A collision will entangle the states of the tracer particle with the buffer gas particle. 
However, we are interested in describing the ensemble of tracer particles only. This is 
sufficient for the description of most experiments which do not probe entanglement between 
the tracer particles and the bath. It is possible that the tracer particle may encounter the 
same bath particle more than once but the bath particle will have undergone further collisions 
entangling its state with those of other bath particles. There is no mechanism which can 
systematically maintain coherence between the tracer particle and the gas particle state 
between successive encounters of the same particles. Thus, for each collision we may treat 
the two-particle state of the colliding pair as the separable state p pair = p <g> p gas , and the 
tracer particle state after the collision is given in terms of the incoming state p by 

p' = tr gas [S(p®p^ s )Si] (9) 
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where S is the two-particle scattering operator which acts on the space l-i v <8> "H rc i <8> "Hcm, 
and tr gas indicates the trace over the buffer gas particle degrees of freedom. Collisions are 
classified as either elastic, in which the internal state does not change, or inelastic, in which 
it does. In the absence of trap loss, the total number of molecules is conserved. Hence, 
probability is conserved and the S operator is unitary. 

The one-particle scattering operator Sq acts on the space 1-L V ® T-i r e\ and 35 1 



S = 1 C M <8> 5*0 . 
The two-particle operator T is related to S by 

S = l + iT. 

Similarly, 

So = 1 + iTo 

where T = 1cm ® Tq and Tq is a one-particle operator with matrix elements 

777 

(p'v'\T \pv) = -^5(E'-E)f u , v (p',p) 



(10) 



(11) 



(12) 



(13) 



where f u 'v{p' , p) is the scattering amplitude for the transition (p', v') ^— (p, v) and = 
mMjim + M) is the reduced mass of the colliding pair. E' = p' 2 /(2m*) + e v > and E = 
p 2 /(2m*) + e v are the total energies for the states (p', u') and (p, v) respectively. 

The time evolution of p for tracer particles in free space undergoing collisional decoherence 
with buffer gas atoms is 

dp\ i 



dt 



coll 



h 



[H™ + H n ,p]+£p 



(14) 



where (dp/dt) co \\ arises from collisions and £ is a dissipative Lindblad operator 
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22 



23] . For a tracer particle with internal states, Vacchini and Hornberger present a heuristic 



derivation of Eq. (1T41 finding [23] 



vv 



/ d 3 p p(p)Re \f vv , (rel(p,P),rel(p,P))j ® \u) (15) 



and 



Ac 



d 2 k 



e iQ A/n L(k, P; Q, Ae)ptf(k, P; Q, Ae) e - lQ ^ 



in 



pL\k, P; Q, Ae)L(k, P; Q, Ae) + tf(k, P; Q, Ae)L f (k, P; Q, Ae)p (16) 
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where 



L(p,P;Q,Ae) 



^2 fw' ( rel (P±Q, P ±q) - y + Q, rel(p± Q , P±q) + ^ + Q 



vv 

Ae=e u -ei 



Q 2 /m* 



m Q 



m . 



A* P±Q + — T + 77 P IIQ + 



Ae 



Q 



2 ' Q 2 /m* 

,1/2 



<8> 1 1/) . (17) 



m» 2 M~ " w Q 2 /m* 

In Eqs. f|T5|) to (|17|) . P and R are respectively the momentum and position operators for 
the tracer molecule, 

is the relative momentum between a buffer gas atom and a tracer molecule with respective 
momenta p and P, and Phq and P_lq represent the components of P respectively parallel 
to and perpendicular to Q (with similar notation for P±q). The integration with respect to 
k in Eq. (TL6]) is over the two-dimensional space perpendicular to Q. 



B. Trap loss in inelastic collisions 

We now include trap loss due to inelastic collisions. The S operator may be decomposed 
into parts corresponding to elastic and inelastic collisions, as may p' as given in Eq. ([9]). We 
define the elastic part of the S operator 

V 

where S uu ' = (v\S\v') is an operator on % re \ <g> %cu- The inelastic part is 

S' m = S- S el = • ( 20 ) 

v^v' 

Unlike S, the operator S el is not unitary. 

The vv' element of p(t) is p uu >(t) = (v\p(t)\v'), which is an operator on the tracer particle 
translational state space. The elastic term in p' is 

/) /el = tr gas [S cl (p®p gas )5 elt ] (21) 

and the vv' matrix element is 

p':l = tr gas SW[/W ® P gas )Sl v , (22) 
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since Sf v , = S uu /. The inelastic term is 

pin = p_ ~,el _ (23) 

pel correS p 0nc [ s to the portion of the ensemble which remains in the trap after the collision, 
while p' m describes particles ejected from the trap. 
The change of p in a single collision is 

Ap = p'-p (24) 

and we may perform a similar decomposition of Ap into terms due to elastic collisions, 

(Apf = p' el - p (25) 

and inelastic collisions, 

( Ap) in = Ap - ( Ap) el . (26) 

If we now consider the continuous-time evolution of the ensemble of trapped tracer 
molecules, the number of molecules in the trap will be continually decreasing due to trap 
loss. Experimental measurements can only probe the molecules in the trap, hence we are 
interested in modelling only the portion of the original ensemble that remains in the trap. 
The unnormalized density operator p el (t) describes the collection of molecules in the trap at 
time t. Since all of the sample is trapped at t = 0, p el (0) = p(0). The fraction of the initial 
sample still trapped at t is trp el (i), where the trace is over the internal and translational de- 
grees of freedom. The molecules ejected from the trap cannot be realistically modelled, nor 
would it be of interest to do so since those molecules cannot be probed by experiment. The 
ejected molecules vanish from the continually depleting trapped ensemble described by p el 
but beyond this their state is not defined. The density operator p(t) , which would ordinarily 
describe the state at t of the entire ensemble of molecules, is therefore not well-defined and 
is not a quantity of physical interest since the outcome of all measurements on molecules in 
the trap is described by p el {t). 

The reduced density matrix describing the internal state of the trapped ensemble is 
Hermitian and has matrix elements pf u ,(t) = Q(2nh)~ 3 J d 3 P p^ u , (P, P; t). These are com- 
plex numbers and should not be confused with the operators p^ v ,(t). With this definition 
/C(0)=/W(0) (cf. Eq. ©). 
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C. Measures of coherence 



There are two quantities of interest for characterizing the degree of coherence between 
the internal states v and v' in the trapped ensemble. The first is \pf u ,\ and the second is 

(el el \ V 2 
(27) 
PvvPv'v 1 J 

These are both nonnegative real numbers, and < i] uu , < 1. For trapped molecules in a 



pure state \pf u ,\ = a/ PwPu'v' anc ^ Vuu' = 1- When there is a complete absence of coherence 
between the states v and v' both \pf u ,\ and i] uu , are 0. Consider an experiment performed on 
the molecules in the trap which measures an observable A = {A UI/ /} that involves only the 
internal but not the translational degrees of freedom. The terms in the expectation value (A) 
that depend on the coherence between the states v and v' are (p^ v ,A v i v + pf, u A ul/ ,)/tr p el = 
2\pf v , \\A vv i | cos 0/tr p cl for some phase angle 0. Coherence in an experiment is often observed 
as temporal or spatial oscillations (i.e. interference fringes) in the value of some observable. 
Hence, r\ vv i represents the ratio of the amplitude of observed interference fringes between v 
and v' to their maximum possible amplitude, which occurs when there is perfect coherence. 
\pf v , | is relevant for applications such as quantum computation where the size of the ensemble 
is an important aspect of the coherent signal. Our goal is to determine \pf u /{t)\ and r] uu i(t), 
and we proceed by first considering the time evolution of the density operator pf u ,(t), which 
describes translational as well as internal degrees of freedom. 

D. Master equation 

In a collision with incoming density matrix p cl , the change Apf u , = p"^ v , — p^ u , is 
Ap^ = tr gas l l -[(T uv + TlWl ® P gas ) - {pt, ® P gas )(^v + Tj v )] 



1 

2 



T w {ft> ® P SaS )^ v , , (28) 



. v" 

where we have used the relation 

l (T uv -Tl) = -Y J Tl»Jv"^ (29) 

v" 

which follows from the unitarity of the S matrix. 
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In the momentum representation Eq. (1281) takes the form 

A P : 1 i/ (P,P') = ^^{p c ^(P,P') J d 3 pp(p) ^((rel(p ! P)HT + T t |rel(p,P)z/) 

- (rel(p, P>' |T + rj|rel(p, 

-\j rf 3 QE(^ rel (p-Q' p + Q) z/ "i T °i rel (P' p »i 2 

+ |(rel(p - Q,P' + QK|T |rel(p,P'K> 
+ J d 3 pp(p) J rf 3 Qpt (P-Q,P'-Q) (rel(p-Q,P)z/|T |rel(p,P-Q)z/) 

x(rel(p,P / -Q)z/ / |r o t |rel(p-Q,P / )z/')} • (30) 

Eq. (130]) contains terms of the form 

(pHTo + Tjipv) = ^5(0)^(p, p) (31) 



and 



(pV'lTolpz.) | 2 = j^ 2 5\E' - E)\fMv',v)? (32) 



which contain the undefined quantities 5(0) and 5 2 (E' — E). Hornberger and Sipe [21] present 
a normalization rule applicable to the single-collision expression Eq. (13"0"|) . However, we use an 
ultimately equivalent technique proposed by Adler in which the normalization is carried 



out simultaneously with the passage from the single collision expression to a continuous-time 
at i 



equation for ^-p el (P, P'; t). 



According to Adler's method, on the left-hand side of Eq. (1301) we replace Ap„ u ,(P, P') 
with Ap^,(P, P') = p^,(P, P'; t + At) - pf„,(P, P'; t), which is the change in the density 
matrix due to collisions during a finite time interval At with a bath ensemble of one particle. 
The coarse-graining time At must be longer than the duration of a collision but short 

n n 

compared to decoherence timescales [211 . l24j . The right-hand side is multiplied by N, the 
number of bath particles in the system volume fl. 

Adler shows that the square of the delta function can be written as 5 2 (E' — 
E) = S(E' - E)5(0) and that the appropriate choice for 6(0) is 5(0) = At/(2nh). 
Hence 5 2 (E' — E) — 5(p' - p)At/(4ir 2 h 2 p). After dividing by At one equates 
[pH(P,P';t + At) -p^,(P,P';t)] /At = [fp^(P, P'; t)] coll , the collisional contribution 
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to 



h \ 2M 



5/2 



2M 



e u ,)p? 1/ (P,P';t) + 



|p c ,l(P,P';t) 



(33) 



coll 



with the other term arising from free Hamiltonian evolution. 

Applying these manipulations to Eq. ( 1301) thus converts, on the left-hand side, 
Ap^,(P,P') — > ^pf u ,(P, P'; £)] . On the the right-hand side the effected conversions 
are 

(2vr/i) 3 



N- 



pu\T + Tjlpi/) -> 47rn^Re/ w (p,p 



iV 



(27rfi) 3 



|(pV'|T |pz,)| 2 ^5 fp' - J p* + 2m*(e v - €„,>)) | Wp', P ) 



(34) 
(35) 



where n gas = AT/fi. 



This leaves the final term on the right-hand side of Eq. ( 1301) . which contains the expres- 
sion (rel(p - Q, P)z/|T |rel(p, P - Q)u) (rel(p, P' - Q)z/|T f |rel(p - Q, P') v') in which the 
initial momenta are different for the two matrix elements. For the decoherence of internal 
state superpositions we are interested in the reduced density matrix for internal states having 
matrix elements p^ v ,(t) = J d 3 P g vv <(P,t) where g vv >{P,t) = (27ih)- 3 Qpf u/ (P,P;t). Eq. (130]) 
established that evolution of the elements of the main diagonal g ul/ >(P;t) is independent of 
the off-diagonal elements pf u ,(P,P';t), P ^ P'. It is therefore sufficient to consider the 
restriction of Eq. (130]) to the P = P' case, thereby avoiding concerns about the extension 
of Eq. (I35p to off-diagonal terms. With straightforward changes of variable, the resulting 
equation for full evolution is 



3 % 
—g vv ,(P,t) = ~{e v , - c 



)g vv ,(P, t) + ^ (1 + r) 3 ( W (P, t) [ d 3 pp((l + r)p + rP) 

{ J 

x 27r*Re(/ w (p, p) - /„v(p, P)) - f (^(P) + ^(P)) 
+ /^ (r(P+pA) + p ^^ 



(36) 



where we have defined r = m/M. The total cross section for scattering in state v with 
pairwise relative momentum p, 

a/p 2 + 2m*{e v - e u »] 



V 



fv»v ( Vp 2 + 2m *(^ - e^")n, 



(37) 
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is related to the scattering amplitude in the forward direction by the optical theorem 
cr* ot (p) = (47rh/p)lm f vu (p, p), which has been used in Eq. fl36|) . 

The evolution of p o1 is not trace-preserving because of trap loss in inelastic collisions. 
Hence there is no equation in Lindblad form for p el such that Eq. fl36l) can be obtained as 
its (Pis, Pis') matrix element. 

With s-wave scattering only, f uv (p',p) — fw(v\ At low momenta the s-wave scattering 
amplitude may be expanded as 

fw(v) = ~ a v + Kp + c u p 2 + ... = -(a v - V ) + (K + i>K)P + ■■■ , (38) 

where the coefficients are in general complex and a v is the complex s-wave scattering length 
for particles in state v j36|. The sign of a v and the notation a v and — f3 u for the real and 
imaginary parts of the scattering length are conventional. 

The elastic cross section is cr^(p) = J d 2 h \f uu (ph, p)| 2 and the inelastic cross section 
°t n (p) = a l ot (p) — ^(p)- At leading orders in p, with the assumption of s-wave scattering 
only, 

A'rrh 

°T(v) = —<A + Kp + d v p 2 + ...) (39) 
af(p) = 4tt [\a u \ 2 + (-a v b* v - a* u b u )p + . . .] (40) 
C(p) = 4tt (h^ + W v - \a u \ 2 + . . j . (41) 

Because a cross section cannot be negative, (3 U > 0. 
We will also use the expansion 

\fw(p) ~ fvv'(p)\ 2 = \ov - <v| 2 - 2Re [(a v - a v ,)(b* u - 6*,)]p 

+ {\K~ K<\ 2 - 2Re[(a u - a vl )(c* v - c^jjp 2 + .... (42) 



III. TEMPERATURE DEPENDENCE OF DECOHERENCE 

At this point it is convenient to introduce dimensionless variables. We introduce a char- 
acteristic length I, which we leave unspecified, that will drop out before the final results. 
We define a dimensionless temperature 9 = 2mkBTl 2 /h 2 , Q = Pl/(Hy9), r = thn gas l/m 
and 7„j/(Q,t) = (h\^9/l) 3 g uu i (^h^/9/l,T[hn gas l/m]^ lS ^. These scalings have been chosen 
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so that the independent variables describing the particle masses are m and r = m/M, a 
property that will be used in Sec. IIV1 
In the scaled variables Eq. (136]) is 



^:7^'( r ) = Ghuu'(r)] 



(43) 



where we have suppressed the Q-argument thereby indicating that we are considering the 
function 7(7") = 7(-, r) as an object which may be distinguished from its value evaluated at 
Q, 7(Q, r). The linear operator G acts on a function /i(Q) as 

G[h](Q)= l { \;^] m h(Q) + (l + rr 



x I HQ) / d 3 q2ni 



f 1 w 1/2 ~ fie 1 / 2 „\ f* ( fie 1 / 2 „ fie 1 / 2 „ 

■/f" I — r~ 1' — r~ q j vi/' I ~ r~ q ' ~~ q 



6 1 ' 2 / d 3 g( i 2 n 



7i(Q - q + gn)g 



; -[rQ+(l+r)q] 2 
^372 

e -(rQ+ri}n+q) 2 



7T 



3/2 



(44) 



G depends on z/ and z/ but we suppress this in the notation for clarity. The initial conditions 
implied by Eq. Q are 

7(Q,0) = p^(0)Q 3/2 e" rQ2 . (45) 
Applying the assumption that there is only s-wave scattering, we substitute Eq. (1381 into 
Eq. (144j) . Since p = Hy/dq/l, we obtain an expansion in powers of 6 1 / 2 : 

G = G + Q XI2 G X + 6G2 + ... , 



where 



with 



vv • 



G [h)(Q) = K^'h(Q) 
{e v — e v ')m 2iti{a v - a*,)(l + r) 



(46) 



(47) 



(48) 



and 

G 1 [h}(Q) = (l + r)U2n 



, . ifr(K-b:,). f ^ f, 3 e -^ 1+r W 



MQ) / d'qq 



71 



3/2 



(i 11 ' I „-(rQ+rgn+q) 2 

+ ^ / fqiP&qKQ-q+tpi)?—;^ } , (49) 



P 



7T 
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and 

.h 2 ( Cu -C*A f „ oe -[rQ + (l+r)q] ; 



G 2 [/i](Q) = (l + r) j fe(Q)2n Mc - 3 C:>) J d 3 qq 2 - 



7T 



3/2 



fj(—n h* — h n* ) f „-(rQ+r(jn+q) 2 

+ ( % V v ) J d\ d 2 h h(Q - q + 9 n) 9 2 ^ > . (50) 

For notational brevity we define J h = J d 3 Q h(Q) for any function h(Q). The zeroth- 
order truncation of Eq. (H3|) has solution 7i, i /(t) = 7 w '(0)e K »" T and yields p ul/ '(r) = 
P«/(0)e K o" T . Observe that 

Re/^' = -2tt(^ + /3,0A <0, (51) 

hence e K °" r is a decaying term except when equality holds, which occurs when the leading 
order of inelastic scattering vanishes for both internal states. The solution of Eq. (|43p has 
the form 7 w /(r) = e K °" t F uu >(t), where 

^F wI {t) = e^GtlF^) + ^G 2 [F^] + . . . . (52) 
ar 

Fw' if) in general does not describe simple exponential decay. Even if we truncate the right- 
hand side at order # 1//2 , Eq. (152j) is difficult to solve because of the complicated form of G\. 
Instead, we adopt a perturbation approach to the analysis of Eq. (H3|) . expanding 

7^ = if + 1/2 li u ' + Olf + ■■■■ (53) 



Substituting this and Eq. (H6|) into Eq. (H3l) . we obtain 

ILtf + O 1 ' 2 ^' + ... = Gohf] + O^iGohf] + Cibo 1 '']) + • • • • (54) 

The initial conditions are 7q i/ '(0) = 7i/^'(0) and 7^ (0) = for k > 1. Equating coefficients 
of like powers of O 1 ^ 2 and solving at the lowest three orders we obtain 

7o -'(r) = e^V,(0) (55) 
7r '( r ) = re ^'^ l[7w;(0 )] (56) 

lf(r) = e K «' T (tG 2 [ 1vv ,(0)) + i^GiIGi^CO)]]) • (57) 

These are the leading terms in the expansion of e Gr p)w(0)], which is the formal solution to 
Eq. (H3J). 
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We substitute Eqs. (155]) and (156]) into Eq. (153|) . and integrate over Q. The integrals in 
the resulting expression are 

_ 2(1 + r) 1 / 2 (2mh{b v - &*,) + 4™^*/ 



Co 



G 2 [7, 
Gi [Gi[ 7 o w '(0)]] 



/o 



7rV2 V. I 2 

3 [2rf( Cl , - c*,) - ^h{a v bl' + 



1 

7T 



n 2 



X 



2/ 3 

3(2r + l)^ + 1 + 2r + 3ra Bin-^ r 



/w(°) 



r r + 1 

The evaluation of Eq. ( |60|) is discussed in the Appendix. 

Converting back to unsealed variables and using Eqs. (I58|) -f l60|) . we obtain 



1 + T 1/2 zft + T ( z£t + 2%£- 



where 



,/ _ i(e u > - e v ) 2iri(a u - a* u ,)hn gas 



B 1 l gas 



m 



1/2 



ih{b v — b* v ,) + 2a v a* v \ 



Z 2A = Qrcn gas k B r 3/2 [ih(c u - c*,) - 2(a u b* v , + b u a* u ,)] 



and 



gas 



"2,2 



m 



3(2r + 1) 1/2 + 



1 + 2r + 3r 2 . _ x r 



sin 



r + 1 



Note that 



2irhn 



gas 



(58) 
(59) 



^(0). (60) 



(61) 

(62) 

(63) 
(64) 



- + 2a,<,] 2 . (65) 



(66) 



Referring to Eq. ( |4~TT) . we observe that equality occurs when the coefficient of the p 1 term 
vanishes in both a l ^(p) and a^(p). As well, 



o5/2 1/21.1/2 



1/2 



Qj/ — cv 



2 )- 



(67) 



Note from Eq. (l4ip that hb l u — \a u \ 2 is the coefficient of the p-independent term of cr™(p). If 
these coefficients are nonegative for erjf and cr™ then Re zj'" < 0. The nonnegativity of a m 
is not sufficient to establish the nonnegativity of the p-independent coefficient. On physical 
grounds the inequality is expected to hold. 
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A. Total coherent signal 



The first of the two measures of coherence between the states \u) and \v') we calculate is 

(68) 



\pUt)\ = ^/p?At)p e X(t). 

From Eq. (l6Tj) and the small- 1 a; | expansion (1 + x) 1 ^ 2 = 1 + x/2 — x 2 /8 + ... we find 



\p?At)\ = \p»Ao)\e-^ 1 



1 + cr ' T i/2 t + T ( c ^ t + c »A_y 



with 



Co 


= -Re*™ 


Ci 


= Re zf 


Ca,i 


= Re z^ 



and 



C 2 , 2 = Re^ 2 ' + (lmzr')' 



(69) 

(70) 
(71) 
(72) 

(73) 



Differentiating Eq. ( 1691) gives 
J t \pUt)\ = \Puw(0)\e- Cot 



x 



-Co + G(l-Cot)T 1/2 + T 



C2,l(l-Co*) + C2,2 (t 



Cot' 



+ ... 



(74) 



The T-independent leading term is the result of trap loss from inelastic collisions. 

It is interesting to determine when the T^-dependent term may be neglected in com- 
parison with the T-independent term. This occurs when 



|Ci(l-Cot)|T 1 / 2 «Co. 
This condition is met for a time interval beginning at t = when 

1/2 „ _Cp_ 



ICi 



and remains valid while 



t < 



| Cl | T V2 • 



(75) 



(76) 



(77) 
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0.5 

T 1 o 

0.25 





12 3 

So* 

FIG. 1: The curves are |&||1 - Co*|^ 1/2 = O.ICo (solid), T 1 / 2 = O.lCo/ICil (dashed) and t = 
O.ldCilT 1 / 2 )- 1 (dotted). 

The regions of the (t,T 1//2 ) space satisfying Eqs. ( 175|) -( 177|) are shown in Fig. [TJ 

It is noteworthy that the T 1 / 2 term becomes significant given sufficiently long time rather 
than remaining negligible indefinitely. Eq. (l43j) does not describe simple exponential decay. 
In the perturbation formulation we may attribute the increasing 6* 1//2 term to the cumulative 
effect of the O 1 ^ 2 component G\ of the time evolution operator G. The ^-independent com- 
ponent Go °f the operator G of Eq. (j4"3"|) causes the pseudodistribution 7 I/I /(Q) to be scaled 
by a Q-independent factor, but there is no change in its shape. G\ contains two terms (see 
Eq. (j49p ). the first of which corresponds to Q-dependent decay. The second term describes 
Q-dependent decay and also the "redistribution" of 7 W '(Q) between different Q values due 
to thermalization, indicated by the Q-dependent convolution appearing in the term. These 
effects change the shape of 7 W '(Q) as opposed to scaling it by a constant. Their cumulative 
contribution will be negligible at short times but significant at sufficiently long times. 
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B. Relative coherence in trapped sample 



The time-dependent relative coherence between the states \v) and \v') in the trapped 
sample is found by substituting Eq. (16Tj) into Eq. (12TI) . We use the small- \x\ expansions of 
(l+x) 1 / 2 and (1+x)" 1 / 2 = l-x/2+3a; 2 /8+. . ., and observe that (z% v ' +z%' v -z% v -z%' u ')/2 = 
and pliX0)pl^(0)puu /2 (0)p- 1 J, 2 (0) = ^(0). We obtain 



where the real-valued coefficients 



i + T y^ 1 + T^ 2jl t + t -^+... 



(78) 



<*\ \ <*\ >6\ 6\ 



— la,, — a,, 1 



1/2 



£2,1 



6,2 



z 2,l "T ^2,1 "T z : 



2,1 ^2,1 



f27rn gas A; i jr 3/2 Re[(a J , - a i/ /)(6* - &*,)] 

z 2.2 ' z 2,2 z 2,2 z 2.2 



(79) 



(80) 



+ 



8nkBnl 



3(2r + f) 1/2 + 



f + 2r + 3r 2 



sm 



r + f 



-4(l + r) 



x 



{ - £v)| 2 - 41m [h{b u a u - b u <a v >){cL u - a*,)] + 2^ 



a,,/ a 



^ + <v| 2 } 



+ 



327r 3 n 2 as fcB , 



0„ — Ct,/ 



(81) 



£x < 0, while the signs of £2,1 and £2,2 are unconstrained. Observe the appearance of 
coefficients from Eq. (@2]) in Eqs. ([79]) and f lHUj) . 
The decoherence rate is 



^'(*) = ^'(°) [ Tl/2 6 + T (6,i + #2,2) + •••]• 



(82) 



We wish to determine when the T-dependent contribution to decoherence as measured by 
i] uu i (t) may be neglected by comparison with the T 1 / 2 -dependent contribution. The condition 
that must be satisfied is 

T 1/2 |6,i + ^2, 2 |«|6|, (83) 
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from which it follows that it is sufficient that both 

rpl/2 ^ 16 _ 

" I6,i| 34 /2 ^i /2 r 3 / 2 |Re [{a v - <v)(&* - K>)\ 

and 



t < 



16 




16,2 


TV2 



(85) 



With Eq. holding, Eq. (185]) is equivalent to the condition that t be comparable to or 
smaller than |6,i|/l6,2|- Fig. |2] illustrates the regions corresponding to Eqs. (l8~3l)-(l85l) for the 
case where £2,1 and £2,1 have the same sign. In obtaining Eq. (1851) . unity has been neglected 
m comparison with T 1 / 2 |5 i2 |/|6,i I which results in a less restrictive condition, indicated in 
Fig. |2]by the solid curve lying above the dotted curve. This is permissible as these conditions 
are order of magnitude estimates rather than strict inequalities. 



23/2^3/2, 



a„ — a,, 



0.5 



0.25 



1 1 1 1 

12 3 

^2,2^2/ 

FIG. 2: The curves are T l / 2 \&,i + ^2,2! = 0.1|£i| (solid), T 1 / 2 = 0.1|£i|/|6,i| (dashed), and 
t = 0.1|6|/(|6,2|T 1 / 2 ) (dotted). The case depicted is that in which £2,1 and £2,2 have the same 
sign. 
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IV. LIGHT BUFFER GAS PARTICLE 



Helium is a likely choice for the buffer gas in experiments, so the low-r regime, corre- 
sponding to a light buffer gas particle, is of particular interest. Since there are two inde- 
pendent mass variables, the r — > + limit may be approached along different paths in the 
two-dimensional space. The behavior as r — > + with M held constant is highly singular. 
We consider the limit as r — > + with m constant, which will allow us to obtain low-r ap- 
proximations and conditions for their validity. The reduced variables introduced in Sec. II III 
were chosen to make r and m independent variables. 



Beginning from Eq. 



we make the low-r expansion 



G = t {€u €u '] m +G + rG 1 + r 2 G 2 + 



(86) 



where for a function h(Q) 



G k [h](Q) = h(Q) d 3 q 



+ 6 1 ' 2 / d 3 qd 2 n 



2m 



ivy \ i H) i " 



f * ( ne 1 / 2 r . he 1 / 2 

Jv'v' \ I H? i M. 



P 



7T 



3/2 



<? fc (Q,q, q) 



g^£fc(Q, q, n)h(Q - q + gn) 



(87) 



and 



So(Q,q,n) = l (88) 
&(Q, q, n) = 4 - 2q ■ (Q + gn) (89) 
&(Q, q, n) = 6 + [2q • (Q + gn)] 2 - (Q + qn) 2 - 8q ■ (Q + gn) . (90) 

Due to the complicated structure of the operators G k a perturbation expansion in small r 
does not lead to anything useful. 
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Inserting Eq. (186]) into Eq. (H3j) and integrating over Q we obtain 



d 

dr' 



{e v ,-e v} m d 

1 6 ' 6 gas fc 



+ / d 3 g 7w '(Q,r) / rf 3 g-^^ (Q,q,q)+^i(Q,q,q)+r 2 6; 2 (Q,q,q) + ...] 



2ttz 



q,^q 



i 



+Q6 1 / 2 J d 2 h 



Suu (ggn, ^q) (^V, 



(91) 



If the quantity in braces is independent of q then q • Q averages to zero over one of the 
angular integrations and the integrals of the terms containing Q in {?fc(Q, q, q) vanish for 
k = and 1. Hence the q integral will be independent of Q. This spherical symmetry 
condition is satisfied at ultracold temperatures. In this case 



d 

dr 
where 



e v > - e u )m 



PtAr) + rwfpfAr) + r 2 I G 2 [^(r)} + . . . (92) 



2tu 



f ( M 1/2 ~ _ f * ( he 1 ' 2 „ fig 1 / 2 „ 

— q>— — qj Ju'u' i - — q> - — q 



^qd 1 ' 2 j d 2 h 



I 



I 2 



(93) 



The evolution of pf u /(r) is exponential up to the first order in r. 

At ultracold temperatures we may use Eq. (138]) . We obtain the expansions 

g i/2 l7r i/ 2 W^-^) + 2a ^] 

qo^ \ %h \ c v ~ c 2> ~ ZKaybl, + b v al,)} | 



i\a v a , 
- 271 — 



-2ir 



t{a u - a u , 



63n 



\ifi 2 {c v -c* u ,)- 2h(a v b* v , + 6„a*')] 
I 3 



+ 



(94) 
(95) 



A. Total coherent signal 



Using the relationship 

i\ P Ur)\ 



2|/C(r)| 



p C X(T)j t pt,(r) + p^ v) {r)— P B X 



T 



(96) 
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and Eq. (192]) we obtain 

^\pt>(r)\ = Rewf\pU + rRewf\ P U + j^Re (fi y J G 2 [ 7 ^)]) + • • • • (97) 

The complicated form of the operators Gk prevents us from finding conditions for the validity 
of the first-order approximation. We will establish conditions for the validity of the zeroth- 
order approximation. 

The r-dependent term is negligible relative to the r-independent term when 



r < 



Rew^ 



Rewf 
1 + (k B T) 



1/2 2 3 / 2 m 1 / 2 [^(&t + tt v ,) - 2(cw + AAQ] 
6m[n 2 (4 + 4,) + 2hRe(a v bt> + 



+k B T 



+ 



(98) 



The coefficient of the T 1//2 dependent term in Eq. (198]) is positive (except in the special 
case when both the numerator and denominator vanish, which we do not consider further), 
thus at low temperatures the mass ratio r at which the zeroth-order approximation is valid 
increases with increasing temperature. The T 1 / 2 -dependent term is insignificant with respect 
to the T-independent term when 



2^m^[h(bl + V vl ) - 2(cw + duM\ ' 



(99) 



When Eq. (199]) holds then the condition for validity of the zeroth-order approximation, 
Eq. (f9"g|) becomes 



r < 1 . 



(100) 



The exponential decay constant for the zeroth-order approximation for the unsealed time 
variable t is then 



Ai 



flTlgagl vv i 

— Rew ~ n gas 



(101) 



m m 

where the validity of Eq. (I9"9"]) has enabled us to drop the T 1,/2 -dependendent and higher 
terms in the small-T expansion of Ai< The factor hn gas l/m comes from the conversion from 
the t time scale to t. At zeroth-order in r, the decrease in \p„ u >\ is determined entirely by 
the leading order terms of the inelastic scattering cross-sections. 
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It is possible that, as in Sec. HTT1 there is a cutoff time for the validity of the zeroth- 
order approximation because of the non-exponential forms of Eqs. (192|) and (}9Tj) . However, 
because the complicated forms of the Gk preclude solving for 7 I/I /(i), we cannot estimate the 
length of time for which the approximation is valid. 



B. Relative coherence in trapped sample 



Differentiating Eq. (j27j) and using Eq. ([92]) . we obtain 



d 



+ r 



2 I / G 2 [ 1w ,{t)\ , / G 2 [ 1v , v {t)\ J G 2 [ lw ]{r) J G 2 [ 1v , v ,]{t) 



Pvu' (t) 



Pv'vir) 



Pvvij) 



Pu'u'ij) 



t)w> (r) + . . . 
(102) 



where 



" " - wl v - w v ' v ' 



e 1 ' 2 



2vr 3 /2/2 



d 6 qq{A-2q 2 ) k e- q2 I d z n 



no 1 ' 2 „ he 1 ' 2 

-qn, 



I 



I 



q 



< 



(103) 



for k = and 1. As before, the decay of coherence is exponential up to first order in r 
regardless of the temperature, and we cannot determine conditions for the validity of the 
first-order approximation. We find 



u 



■loj, — H ^— Re [(a„ - a u/ )(b u - b u ,)\ + 



I 2 
GirftO 



(104) 



Z 3 



Re[(a I/ - a v :){b* v - 6*,)] 



K - K,\ 2 - 2Re(a, - 0l/ /)(c* - c*,)]} + • • • • (105) 



The r-dependent term is negligible compared to the r-independent term when 
u 



r < 



-2 djy (2;/ 



37r 1 /2(2mA ;iJ T) 1 /2Re [( a „ - aj/ )(&* - K,)) 

_ 8\a u - a v >\ 2 {\b v - b v >\ 2 - 2Re[(a u - a„/)(c* - c*,)]} 
9tt {Re [(a, - ^ )(&*-&*,)] } 2 



(106) 



25 



The last expression is an expansion at small T 1 / 2 using Eq. (|38|) . which has been obtained 
with the ultracold temperature assumption that there is only s-wave scattering. The terms 
after the first may be ignored when 

2\a v -a v i\ 2 

{2mk B T) 1/2 < 



3wVa|Re[(a„-(v)(65-&^)] 




8\a v — <v| 2 { |6„ — b u ,\ 2 — 2Re[(a„— a v 


')( C J- C *,)]} 


97r{Re[(a„-a„,)(b;-b*,)]} 2 



(107) 



in which case the condition Eq. (1106 j) for the validity of the zeroth-order in r approximation 
becomes 



m 3/2 < 



2\a v - a u A 2 M 



(108) 



^{2k B T)^\Ke{{a v - <v)(&* - K,))\ ' 
The exponential decay constant for in the zeroth-order approximation for the unsealed 



time variable t is 

A 2 : 



w 



n 



rn 



{2mk B T) 1/2 A7T 1/2 \a u - a v ,\ 2 - (2mfc iJ T)67rRe [(a^ - cv)(6* - 6*,)] 

+ (2m/c B T) 3/2 87r 1/2 {\b v - b u ,\ 2 - 2Re[(a u - a v ,){cl - <,)]} + •••)• (109) 



Satisfying Eq. f)107p is not sufficient to determine the order of (2m/csT) 1//2 at which this 
expansion may be terminated. If the denominator on the right-hand side of Eq. fllQ7p is 
comparable to or greater than unity then all terms after the first may be ignored. 

The remarks at the end of Sec. IIV Al regarding possible cutoff times for the validity of 
the zeroth-order approximation also apply to the approximation in this section. 



V. CONCLUSION 



The decoherence of trapped particles prepared in a coherent superposition of internal 
states as a result of collisions with buffer gas atoms at very low temperatures is described 
by Eqs. (171|) and (182|) . The time evolution of the coherence between the internal states v 
and v' is parametrized by the complex scattering lengths of the particles in these two states 
and the other coefficients in the expansions of the s-wave scattering amplitudes f uu (p) and 
fv'v'ip) given in Eq. (1381) . Given these scattering parameters, Eqs. (17^|) and (182]) can be used 
to calculate the decoherence rates taking into account elastic collisions as well as inelastic 
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contributions. These equations and their time-integrated forms, Eqs. (169|) and ( ITS]) , show 
that coherence between the internal states of the trapped tracer molecules, and its time 
evolution, can be represented in the limit of low temperatures as an expansion in powers of 
T 1 / 2 . Regardless of the temperature, in both Eqs. (|74|) and (I52"j) . terms with higher-order 
T 1/,2 -dependence become significant given sufficient time. 

Eq. (|7"4"1) describes the time evolution of the magnitude of the off-diagonal reduced matrix 
element \pf u /\ which characterizes the total coherent signal between the internal states v 
and v' . This quantity depends on the size of the remaining trapped population relative 
to the original number of trapped molecules at t = 0. Eq. (17"4"|) is a low-T expansion of 
JjI/^/I in powers of T 1 / 2 . The leading term in the expansion of Jj|p^/| in powers of T 1 / 2 
is a temperature-independent term arising from trap loss caused by inelastic collisions. It 
is determined by the complex parts of the scattering lengths of the states v and z/, which 
are measures of the low-temperature inelastic scattering cross-section (cf. Eq. ( HTj) ). If 
Pu — 0v' — 0) then the inelastic cross- sect ions vanish at the leading order in the collision 
momentum p and the lowest order term in jjj.\pf u ,\ varies as T 1 / 2 . The T-independent term 
varies as m^ 1 and the T^-independent term varies as mj' 1 hence loss of coherent signal is 
slower for heavier tracer molecules or buffer gas atoms. This is because the mean velocity 
of a heavier particle is slower for a given temperature and hence the collision rate is lower. 

Eqs. f!76p and (I77p give conditions on temperature and time for which the decay of \p„ v ,\ 
may be approximated by the temperature-independent term of Eq. ( 1711) . The regions of 
(t, T 1//2 ) space satisfying these conditions are indicated in Fig. [TJ 

Eq. (152]) gives the time-dependence of i] uu r(t), which is the fraction of coherence between 
states v and v' in the trapped ensemble relative to that of a pure state. The quantity 
7] vv i depends only on the state of the trapped ensemble and not on the size of the trapped 
sample relative to the original population. As a consequence, the expansion of ^kT\ V v' in 
powers of T 1//2 has no T-independent term. The leading order term has T 1 / 2 -dependence. 
It is determined by the square magnitude of the difference in complex scattering lengths 
between the states v and z/, \a u — a u i\ 2 = (a u — a u /) 2 + (/?„ — (3 u r) 2 and therefore depends on 
both elastic and inelastic scattering properties. 

If both of the scattering lengths are equal in their real and imaginary parts, Eq. (I52"j) 
shows that the T 1 / 2 -dependent term of -^i] uu i vanishes. In this case the leading term will 
depend on a higher power of T 1 / 2 . These higher order terms will in general be nonzero 
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since they are governed by the coefficients of terms of order p and higher (e.g. b u , b u >) in 
the expansion given by Eq. (|38|) (see Eq. (HDj) ). The decay of r\ uv i may be approximated by 
the leading order, T 1//2 -dependent term of Eq. (182]) when Eqs. (1841) and (185]) hold. These 
conditions are illustrated in Fig. [2J The T 1//2 -dependent term varies as m 



1/2 



For low bath particle/tracer particle mass ratios m/M, the decay of coherence measured 
by either \p„ v >\ or r\ vv i is found to be exponential up to first order in m/M, regardless of 
temperature. When Eqs. (199]) and fl 1 00 j) are satisfied, the zeroth-order approximation for the 
evolution of \p^ u /\ may be applied. In this approximation decay is exponential with decay 
constant Ai given by Eq. (jlOip . The zeroth-order approximation, with exponential decay 
with decay constant A2 given by Eq. (I109p . may be applied to r\ vv i when Eqs. fll07[) and 
QUID hold. 

The scattering lengths a v for ultracold atoms are usually found from thermalization 



measurements 



37H39], which yield results with significant error bars. Decoherence rates can 



be measured sensitively by observing the damping rate of coherent oscillations |27l . |28| . |40|. 
The imaginary parts of the scattering lengths (3 U are relatively easy to find precisely by 
trap loss measurements. Since £1 and A2 (at low temperatures) are proportional to [a v — 
a u') 2 + ifiv — A/) 2 , Eqs. ( l82l) or Eqs. ( 11021) and (11091) could therefore provide the basis for a 
method of precise determination of the real part of the scattering length a v < when the other 
scattering length a v is known. Precise measurements of scattering lengths in ultracold gases 
may be used as a probe of fundamental constants and symmetries in nature 41] . 
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Appendix A 

Here we evaluate the integral in Eq. (160]) . From Eqs. (l4"5i) and (T49|) . and using 7q i/ (0) = 
7^/(0), we obtain 



G 1 [G 1 [ 1 f (0)}} 



2mh(b v -&*/)+ ^a u a* u , 



} r/w (0) (Al) 



r 3/2 7r 9/2' 
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where 

2 



A(r) = J d 3 Qe-^ d 3 qq e -^ + ^ . (A2) 



Integrating with respect to Q yields 



'M- ^y^ fe.' (A3) 



where 



Ai(s) = J d 3 xd 3 yxyexp [— (a; 2 + y 2 — 2sx • y)] 



71-2 d 2 . , , . . , N 

7S^ (s) (A4) 



with 



A 2 (s)= dx dy e~ {x +y ) (e 2sxy - e~ 2sxy ) . (A5) 

We convert to polar coordinates (i?, 0) and integrate with respect to R. We then make 
the change of variables u = s(l — s 2 ) -1 / 2 cos 26 1 and integrate with respect to u, obtaining 
A 2 (s) = (1 — s 2 ) -1 / 2 sin -1 s. Performing the differentiations of Eq. (1A4[) . we substitute the 
result into Eq. (IA3[) . followed by substitution into Eq. flAll) to obtain Eq. (160|) . 
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